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ABSTRACT 

We discuss the imprints left by a cosmological evolution of the star formation 
rate (SFR) on the evolution of X-ray luminosities Lx of normal galaxies, using 
the scheme proposed by White and Ghosh (1998, WG98), wherein the evolution 
of Lx of a galaxy is driven by the evolution of its X-ray binary population. As 
indicated in WG98, the profile of Lx with redshift can both serve as a diagnostic 
probe of the SFR profile and constrain evolutionary models for X-ray binaries. 
We report here the first calculation of the expected evolution of X-ray luminosities 
of galaxies, updating the WG98 work by using a suite of more recently developed 
SFR profiles that span the currently plausible range. The first Chandra deep 
imaging results on L^-evolution are beginning to probe the SFR profile of bright 
spirals: the early results are consistent with predictions based on current SFR 
models. Using these new SFR profiles, the resolution of the "birthrate problem" 
of low-mass X-ray binaries (LMXBs) and recycled, millisecond pulsars (WG98) in 
terms of an evolving global SFR is more complete. We discuss the possible impact 
of the variations in the SFR profile of individual galaxies and galaxy-types. 

Subject headings: binaries: close— stars: formation— stars: evolution— galaxies: 
evolution— X-rays: galaxies— X-rays: stars 

1. INTRODUCTION 

The global star-formation rate (SFR) has undergone strong cosmological evolution: it 
was ~ 10 times its present value at z ~ 1, had a peak value ~ 10-100 times the present one 
in the redshift range z ~ 1.5-3.5, and declined again at high z (Madau et al. 1996; Madau, 
Pozzetti & Dickinson 1998, M98; Blain, Small, Ivison & Kneib 1999, B99a; Blain et al. 1999, 
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B99b, and references therein). Details of the SFR at high redshifts are still somewhat 
uncertain, because much of the star formation at 2 ^ 2; ^ 5 may be dust-obscured and so 
missed by optical surveys, but detected readily through the copious submillimeter emission 
from the dust heated by star formation (Hughes et al. 1998; Barger et al. 1999). 

The X-ray emission of a normal galaxy [i.e., one without an active nucleus) is 
dominated by the integrated emission of the galaxy's X-ray binary population (see, e.g.., 
Fabbiano 1995). In WG98 we discussed the effects of an evolving SFR on the evolution of 
X-ray binary populations of galaxies, and so on that of the X-ray emission from normal 
galaxies. We suggested that X-ray luminosities Lx of normal galaxies could show significant 
evolution (by up to a factor ~ 10), even in the relatively nearby redshift range z ~ 0.5-1.0. 
We also showed that an evolving SFR could resolve the "birthrate problem" involving 
LMXB and their descendant "millisecond" radio pulsars (MRP, see Kulkarni & Narayan 
1988; Lorimer 1995, L95). 

The SFR profile used in WG98 — the only profile available at the time — was based 
on the optical/UV data alone. Over the past three years, there has been considerable 
progress in our understanding of cosmic star-formation history. In addition, very deep 
X-ray imaging with Chandra is beginning to detect normal galaxies in the redshift range 
z ~ 0.5-1.0, so that comparison with observations is becoming possible for the first time. 
In this Letter, we consider quantitatively the key imprints left by SFR evolution on the 
Lx-evolution profiles of normal galaxies, using the best SFR models currently available. 
We briefly discuss the recent results of Brandt et al. (2001, BranOl) from the ultradeep 
Chandra imaging of the Hubble Deep Field North (HDF-N). A detailed calculation of 
the expected X-ray flux distribution of the HDF-N galaxies based on the results of this 
paper is reported by Ptak et al. (2001, PtakOl). We also reconsider the resolution of the 
LMXB-MRP birthrate problem using the new SFR proflles, and flnd that the calculated 
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rates for both (a) whole populations of LMXB and MRP, and, (b) short-period systems, 
are consistent with observation for some SFR profiles suggested recently to account for the 
multiwaveband SFR data. We discuss the relative roles of global SFR profiles on the one 
hand, and the profiles of individual galaxies or galaxy-types on the other. 



2. X-RAY LUMINOSITY EVOLUTION WITH EVOLVING SFR 

The total X-ray output of a normal galaxy can be modeled as the sum of those of its 
HMXB and LMXB as per the WG98 scheme, wherein the evolution of each species "z" is 
described by a timescale (see WG98). To study the effects of the dependence of Tj on 
the binary period and other evolutionary parameters, we run the evolutionary scheme over 
ranges of hkely values of Tj given in the literature. The evolution of the HMXB population 
in response to an evolving star-formation rate SFR(t) is given by 

= a,SFR(t) - (1) 
ot thmxb 

where riHMXB is the number density of HMXBs in the galaxy, and thmxb is the HMXB 

evolution timescale. ah is the rate of formation of HMXBs per unit SFR, given 

approximately by ah = | /binary /prim /sn; where /binary is the fraction of all stars in binaries, 

/prim is that fraction of primordial binaries which has the correct range of stellar masses 

and orbital periods for producing HMXBs (van den Heuvel 1992, vdH92 and the references 

therein), and /g^ ~ 1 is that fraction of massive binaries which survives the first supernova. 

In these calculations, we have adopted a representative value thmxb ~ 5 x 10^ yr according 

to current evolutionary models. In our introductory model here, thmxb includes both 

(a) the time taken (~ 4 — 6 x 10^ yr) by the massive companion of the neutron star to 

evolve from the instant of the neutron-star-producing supernova to the instant when the 

"standard" HMXB phase begins, and, (b) the duration (~ 2.5 x 10^ yr) of this HMXB 

phase (vdH92 and references therein). Since the second timescale is neghgible compared to 
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the first, little error is made by approximating this two-step process by a single step with 
an overall timescale thmxb- 

Two basic methods of LMXB production have been discussed. In the cores of dense 
globular clusters, they can be produced by the tidal capture of a neutron star by a normal 
star. Over the rest of a galaxy, stellar densities are insufficient for tidal capture, and LMXBs 
are produced by the evolution of primordial binaries (see, e.g., Webbink, Rappaport & 
Savonije 1983; Webbink 1992). In this paper, we consider only the latter mechanism. For 
spiral galaxies, at least, this must be the dominant mechanism, since the globular-cluster 
LMXB population in such galaxies only accounts for a relatively small fraction of the total 
X-ray luminosity. 

LMXB evolution from primordial binaries has two distinct stages (WG98) after the 
supernova produces a post-supernova binary (PSNB) containing the neutron star. First, 
the PSNB evolves on a timescale tpsnb due to nuclear evolution of the neutron star's 
low-mass companion and/or decay of binary orbit due to gravitational radiation and 
magnetic braking, until the companion comes into Roche lobe contact and the LMXB turns 
on. Subsequently, the LMXB evolves on a timescale tlmxb- Since tpsnb and tlmxb are 
comparable in general, we must describe the two stages separately (WG98) by: 



Here, npsNB and nLMXB are the respective number densities of PSNB and LMXB in the 
galaxy, and ai is the rate of formation of LMXB per unit SFR, given approximately by 
ai = I /binary /prim/sN> whcrc /binary is the fractiou of all stars in binaries, f^^^^ is that 
fraction of primordial binaries which has the correct range of stellar masses and orbital 
periods for producing LMXBs, and fgj^ is that fraction of such binaries which survives the 
massive star's supernova. 




^1 = a,SFR(t) - "^^^ 



(2) 



^?^LMXb(^) _ ?^PSNB(t) _ nLMXB (t) 

dt TpsNB Tlmxb 



(3) 
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We display evolution in terms of the redshift which is related to the cosmic time t 
by tg = 13{z + where tg is t in units of 10^ yr, and a value of Hq = 50 km Mpc~^ 

has been used^. We consider the suite of current SFR models detailed in Table 1 to cover a 
plausible range, using the parameterization ofB99a,b. Models of the "peak" class have the 
form 



These functional forms are not unique, but useful, since they have a convenient low- 2; limit, 
SFR(2;) cx (1 + zy, where all SFR profiles must agree with the optical/UV data (M98), and 
since the model parameters can be manipulated to mimic a wide range of star-formation 
histories (B99b). Peak-class profiles are useful for describing (a) SFRs determined from 
optical/UV observations, i.e., Madau-type (M98) profiles, called "Peak-M" in Table 1, 
and, (b) more general SFRs with enhanced star formation at high z, a good example of 
which is the "hierarchical" model of B99b, wherein the submillimeter emission is associated 
with galaxy mergers in an hierarchical clustering model of galaxy evolution. Anvil-class 
profiles are useful for describing the results of "monohthic" models. The "Gaussian" model 
(B99a,b)is an attempt at giving a good account of the SFR at both low and high z by 
making a composite of the Peak-G model (see Table 1) and a Gaussian star burst at a high 

^For ease of comparison with WG98 and M98, we continue to use here a Friedman cos- 
mology with go = 1/2; other cosmologies will be considered elsewhere. Other values of the 
Hubble constant lead to a straightforward scaling: for Hq = 70 km s~^ Mpc~^, for example, 
tg ~ 10{z + so that the results remain unchanged if all timescales are shortened by a 

factor of 1.3. 





while those of the "anvil" class have the form 
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redshift Zp, i.e., a component 

SFRGauss(^) = eexp |-M_-^} . (6) 

Based on the IRAS luminosity function, this component is devised to account for the 
high- 2; data, particularly the submillimeter observations (B99a). For its parameters (see 
Table 1), we have used the revised values given in B99b. In all models considered here, no 
galaxies exist for sufficiently large redshifts, ^ > 10. 

Figures 1 and 2 show the prompt evolution of HMXBs and the slow evolution of 
LMXBs, and the evolution of the total X-ray binary population, where the two components 
have been so weighted as to represent the total X-ray emission from the galaxy (the 
weight-ratio is the product of ah/oti and the ratio of the average luminosity of HMXBs 
and LMXBs). The HMXB profile closely follows the SFR profile because thmxb is small 
compared to the SFR evolution timescale. By contrast, the LMXB profile has a significant 
lag behind the SFR profile because tpsnb and tlmxb are comparable to SFR evolution 
timescale: depending on these timescalcs, the LMXB profile peaks at redshifts ~ 1-3 later 
than the HMXB profile (as seen clearly in Figs. 1 and 2), which is a characteristic signature 
of SFR evolution (WG98). 

Figure 1 compares the Lx-evolution corresponding to the (Madau or Peak-M) SFR 
profile and the evolutionary timescales we originally used in WG98, i.e., (a) tpsnb = 1.9 
Gyr, Tlmxb = 0.1 Gyr for the whole LMXB population, and (b) tpsnb — 0.9 Gyr, tlmxb — 
0.5 Gyr for the short-period systems. In Figure 2, we display the Lx-evolution for a range 
of SFR profiles — Peak-M, Hierarchical, Anvil-10, and Gaussian, the evolutionary timescales 
being held fixed at tpsnb = 1-9 Gyr, tlmxb = 1-0 Gyr. Between them, the two figures 
thus explore the effects of (a) varying the evolutionary timescalcs for a fixed SFR profile, 
viz. Peak-M, and (b) varying the profile for a fixed set of evolutionary timescales. For 
sufficiently slow LMXB evolution, the galaxy's X-ray emission is dominated by LMXBs at 
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low redshifts {0 ^ z ^ 1), and by HMXBs at high redshifts. As a result, the total Lx-profile 
is strongly influenced at high redshifts by the SFR proflle (see Fig. 2), and may actually 
show hints of a double-peak structure (WG98) at intermediate redshifts {e.g., the Gaussian 
profile in Fig.2). For sufficiently fast LMXB evolution, on the other hand, LMXBs and 
HMXBs may have comparable contributions to Lx at low redshifts (see Fig.l). 

BranOl estimate that the average X-ray luminosity of the bright spiral galaxies at an 
average redshift z ^ 0.5 used in their stacking analysis is about a factor of 3 higher than 
that in the local Universe {z < 0.01). This observed evolution, Lx(0.5)/Lx(0.0) ~ 3, can 
be compared with our theoretical results in Table 2. The degree of evolution from 2; = to 
z — 0.5-1.0 increases from Madau-type profiles to those with additional star formation at 
high redshifts, the numbers for the Peak-M profile being in best agreement with BranOl. 
For a given profile, the evolutionary factor is smaller for a slower evolution of LMXBs 
(roughly measured in this context by the total time tpsnb + tlmxb), as expected. More 
sophisticated estimates of the expected Lx-distribution of HDF-N galaxies, based on these 
evolutionary scenarios, are described in PtakOl. 



3. BIRTHRATE PROBLEM: FURTHER CONSIDERATIONS 

As discussed in WG98, the evolution of MRPs from LMXBs is described by an 
equation similar to Eq. 3, but involving the MRP number density umrp, and MRP evolution 
timescale tmrp ~ 3 x 10^ — 3 x 10^° yr (Camilo et al. 1994). Our evolutionary scheme 
yields the profiles of LMXB and MRP evolution for a given SFR profile, giving the number 
ratio, Nr = nMRp/riLMXB, and the rate ratio, Rr = ^^mrp /!iLMXB ^-j^g ^RP and LMXB 
populations at the present epoch {z — 0). We showed in WG98 that the Peak-M profile 
yielded i?,. — 1 for the overall MRP and LMXB populations, in agreement with current 
observations (L95). However, for short-period systems (LMXB periods ^ 3 days), this 
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profile yielded a value i?^ — 3, smaller than the current value i?^ ~ 8 estimated from 
observations (L95). On repeating these calculations with other SFR profiles which have 
additional star formations at high redshifts, e.g., the Gaussian and hierarchical profiles, we 
now find reasonable agreement for both whole populations and short-period systems, since 
these give i?^ ~ 6 — 8 for short-period systems, for plausible evolutionary timescales. Thus, 
there appears to be no discrepancy between observation and the idea that SFR evolution 
naturally leads to values of Rr which can be ~ 1 as well as considerably above unity (but 
see §4). 

A new development in SFR research since WG98 has been the study of star-formation 
histories of individual galaxies and various galaxy-types. SFR profiles of individual galaxies, 
ranging from those in the Local Group to those in the HDF at redshifts 0.4 ;^ z 1, have 
been inferred, using a variety of techniques (Glazebrook et al. 1999, Abraham et al. 2000, 
Hernandez et al. 2000). For various galaxy-types, models of Spectro-photometric evolution, 
which use the synthesis code Pegase and are constrained by deep galaxy counts, have been 
developed (Rocca-Volmerange and Fioc 2000), leading to a model SFR profile for each type. 
The birthrate problem was originally formulated entirely in terms of observations in our 
own galaxy, so that, strictly speaking, we should use the SFR profile of our galaxy. As none 
of the above techniques can be used in our own galaxy, it is difficult (but not impossible) to 
determine its SFR profile: calculations with such a profile will be described elsewhere. 

4. DISCUSSION 

We have shown that different global SFR profiles within the currently admissible range 
lead to very different Lx-evolution profiles, so that the latter profiles can be an independent 
X-ray probe of cosmic star-formation history. Our results and a more detailed consideration 
by PtakOl indicate that the early Chandra results (BranOl) are consistent with the Peak-M 
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profile: let us now clarify the underpinnings of this result by applying the considerations 
of global and individual SFR profiles summarized above. BranOl used bright spirals for 
their stacking analysis. Rocca-Volmerange and Fioc (2000) have shown that the model 
SFR profile for such (Sa-Sbc) spirals rises roughly in a Madau fashion from z — to z ^ 1 
(which these authors ascribe to a bias in the original sample used to construct the Madau 
profile towards bright spirals), and thereafter flattens to a roughly constant value ~ 12 
times that at ^ = 0, falling again at ^ ^ 7. We can roughly represent this profile in the 
range < z ^ 7 by an anvil-type profile (see §2), with the parameter z^ax as given in Table 
1, and the parameter p ^ 2.7. For such a profile with the timescales tpsnb = 1-9 Gyr, 
tlmxb = 1-0 Gyr, as in Figure 2, our evolutionary scheme gives Lx(0.5)/Lx(0.0) = 3.3, and 
Lx(1.0)/Lx(0.0) = 5.4, in good agreement with both the BranOl results and the Peak-M 
results given in Table 2. We now see why the Peak-M profile appears to give a good account 
of the BranOl results. In effect, the BranOl analysis may be probing the SFR profile of 
only the bright spirals in HDF-N, and the fact that the Peak-M profile is consistent with 
the BranOl results does not imply that the global SFR necessarily follows the Peak-M 
profile. Thus, there is no basic contradiction between the results of §2 and §3. However, the 
fact that our global analysis of the birthrate problem (§3) gives the best agreement with 
observation for SFR profiles (Gaussian and hierarchical) which have larger SFRs at high 
redshifts than the model Sa-Sbc profile of Rocca-Volmerange and Fioc (which, in turn, has 
a larger SFR at high redshifts than the Peak-M profile) remains to be understood fully. 
The individual-galaxy treatment may hold the key here. 

Signatures of star formation in the relatively nearby {z < 1) Universe have usually 
been sought at wavelengths other than X-rays. In this redshift range, the X-ray probe 
described here is a "fossil record" study of star formation at earher epochs for sufficiently 
slow LMXB evolution (see Fig. 2), since LMXBs then preserve this record for ~ 1 Gyr after 
the star- formation peak is gone. Conversely, for a known SFR profile, this is a unique 
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proving ground for theories of LMXB evolution, which have been built almost wholly on 
the basis of our experience in the current epoch: this new probe is, at present, our only 
direct means of watching LMXB evolution unfold over cosmic time. 

It is a pleasure to thank L. Angelini, R. Griffiths, R. Mushotzky, and A. Ptak for 
stimulating discussions. 
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Table 1: Star Formation Rate (SFR) Profiles' 



Model 


^max 


P 


Comments 


Peak-M 


0.39 


4.6 


Madau profile 


Hierarchical 


0.73 


4.8 


Hierarchical clustering modeP 


Anvil- 10 


1.49 


3.8 


Monolithic models 


Peak-G 


0.63 


3.9 


"Peak" part of composite Gaussian 


Gaussian 


N/A 


N/A 


Gaussian starbursf^ added at high-z 



"Model parameters taken from B99a,b. 



''We have chosen the B99b model with a dust temperature 45 K. Note that the parameter p 
for the SFR equals 3/2 plus the value of the parameter p occurring in equation (16) for the 
merger efficiency in B99b. 

'^Parameters of Gaussian starburst component (see cq. [6] of text), are from the modified 
model given in B99b, i.e., Zp — 1.7, a — 1.0 Gyr, and 6 — 70. 



Table 2: Evolution of X-ray Luminosity Lx 



Model 


TpSNB 


tlmxb 


Lx(0.5) 
ix(O.O) 


Lx(l.O) 
LxiO.O) 


Pcak-M 


1.9 


0.1 


3.9 


5.4 


Peak-M 


0.9 


0.5 


4.6 


6.8 


Peak-M 


1.9 


1.0 


3.4 


4.1 


Hierarchical 


1.9 


1.0 


6.2 


11.3 


Anvil- 10 


1.9 


1.0 


5.4 


10.1 


Gaussian 


1.9 


1.0 


7.5 


16.0 
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Fig. 1. — Evolution of HMXB population (dotted line), LMXB population (dash-dotted line), 
and the total X-ray luminosity Lx (solid hne) of a galaxy with a given SFR (dashed line) . As 
absolute ordinate scales are irrelevant for these evolutionary profiles, they have been adjusted 
for convenience of display: Lx always starts at 0.001 at ^ = 0, so that its evolution can be 
immediately read off the figure, and SFR always starts at 0.01 at 2; = 0, so that different SFR 
profiles can be readily compared. This figure is for the Peak-M profile (see text), showing 
the effect of varying the evolutionary timescales tpsnb and tlmxb, whose values are written 
on each panel. The timescales used here are those used in WG98, corresponding to whole 
LMXB populations and short-period systems. 
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Fig. 2. — Same as Figure 1, but showing tlic effect of varying the SFR profile. The evolu- 
tionary timescales are kept fixed at tpsnb = 1-9 Gyr and tlmxb = 1-0 Gyr for all cases, and 
SFR profiles from Table 1 are used. Each panel is labeled by the name of its SFR profile. 
Evolutionary factors from this and the previous figure are collected in Table 2 and described 
in the text. 



